End-to-end entanglement in Bose-Hubbard chains 
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We study the ground state as well as the dynamics of chains of bosons with local repulsive 
interactions and nearest-neighbour exchange using numerical techniques based on density matrix 
decimation. We explore the development of entanglement between the terminal sites of such chains 
as mechanisms are invoked to concentrate population in these sites. We find that long-range en- 
tanglement in the ground state emerges as a result of transfer taking place across the length of 
the whole chain in systems with appropriate hopping coefficients. Additionally, we find appropriate 
perturbations to increase the entanglement between the end sites above their ground state values. 



I. INTRODUCTION 

As constant experimental developments expand our 
possibilities of materializing a new generation of comput- 
ing devices, the need for a better understanding of quan- 
tum phenomena becomes more important. Outstanding 
emphasis has been placed on the concept of entanglement 
as, for example, a resource for quantum information pro- 
cessing and in general as a tool to study the contrast 
between classical and quantum physics. One fundamen- 
tal problem that is now the subject of intensive research 
is the emergence of quantum correlations between distant 
sites of a quantum system. For most many-body systems 
this is notoriously difficult and entanglement between dis- 
tant sites can only be achieved through clever engineering 
[l[ or non-equilibrium dynamics [1, S (3) S , where, in 
the latter case, either intricate methodology is required 
or the amount of entanglement between the ends, though 
finite, is not substantial in amount. Typically, in quan- 
tum many-body systems, even when total correlations 
are long range, the entanglement between individual con- 
stituents such as spins is extremely short range, such as 
between nearest or next to nearest neighbours Jj. Thus 
it is already interesting when in some many-body system 
an entanglement between its farthest components, such 
as the spins or harmonic oscillators at the very ends of an 
open chain, can be generated and even more interesting 
when it is substantial in amount. Most of the work done 
in this field has been carried out in light-matter systems 




FIG. 1: Potential realization of a Bose-Hubbard chain in an 
optical lattice. The above finite lattice can be realized in 
one cell of a superlattice. The end sites of the lattice are 
depicted differently in the figure as we will require them to 
have vanishing repulsion in our study. 



and spin chains, but more recently, the development on 
new numerical techniques and algorithms [8] has opened 
the possibility of addressing more challen ging scenarios 
such as chains of interacting bosons 0, [^, [HI ■ Usu- 
ally, where bosonic chains have been studied, though true 
quantum correlations have indeed been found to emerge 
between distant individual sites, they are generically not 
"substantial" in magnitude 0, [l^- Thus substantial 
entanglement between distant sites of a Bose-Hubbard 
model, even if with some reasonable engineering of the 
Hamiltonian, should be of great interest, not only be- 
cause of the uncommon nature of such long distance en- 
tanglement, as motivated above, but also because from 
the practical point of view this entanglement is known to 
be "distillable" , i.e., convertible from an impure (mixed 
state) to a pure useful form through local actions. Den- 
sity matrices whose partial transposes display negative 
eigenvalues are distillable [3] ~ so if one finds a mixed 
entangled state with the above property then, one can 
already claim, at least in principle, to have a resource for 
quantum communication. 

Chains of bosons are more often than not well de- 
scribed by the Bose-Hubbard model [H, [H, [11]. In 
this model bosons can hop between neighbouring sites 
while undergoing local repulsion when several bosons oc- 
cupy the same site. The development in atom cooling 
techniques has led to a significant rise in the amount of 
experiments related to the Bose-Hubbard model. The 
transition from a superfluid to a Mott insulator has been 
verified in numerous experiments of cold atoms in opti- 
cal lattices as for example in Refs. [H, \\J, 18, 19, 2(|, 
just to mention a few. As a result, physical realizations 
of quantum systems displaying key features such as long 
range entanglement is becoming more and more feasi- 
ble. Motivated by these advances, in this paper we anal- 
yse from a fundamental point of view the emergence of 
long-range quantum correlations in Bose-Hubbard chains 
(Fig. [T]). There has been one previous investigation on 
long range entanglement in a Bose-Hubbard chain where 
the dynamics is effectively reducible to a single particle 
propagation in a lattice p\. Our study is, however, irre- 
ducibly a many-particle situation and can even result in 
more entanglement than what a single particle hopping 
in a lattice can ever do. It is also pertinent to point out 
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that there has been much interesting work on the various 
forms of entanglement contained in Bose gases [111, , 
but here we concentrate on a specific type of entangle- 
ment, namely that between two individual sites at the 
very ends. If implemented in the laboratory, the entan- 
glement obtained from the scheme presented in this work, 
as it is distillable, would provide a resource for quantum 
information applications. In addition to the end-to-end 
entanglement, which has been motivated above, in this 
paper we will also present some results for the entangle- 
ment between a couple of other bi-partite partitions of 
the chain to fully appreciate the distribution of entangle- 
ment between parts of a finite chain as the strength of 
the local repulsion between bosons is varied. 



II. THE HUBBARD MODEL: 
CHARACTERISTICS AND METHODS 

We assume a chain of size N with M bosons, nearest- 
neighbour hopping and open ended boundary conditions 
governed by a Bose-Hubbard Hamiltonian with variable 
coefficients, 



H = Y^ -Yalak{alak - 1) - X] A(al+iafe + a+^fc+i)- 
fe=i fc=i 

(1) 

Constants Uk and Jk account for the on-site repul- 
sion and hopping respectively while the annihilation and 
creation operators a\. and Ofc obey the usual commut- 
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ing rules [ai , a[] 
Physical realizations of this Hamiltonian include experi- 
ments where remarkable control of the ratio U/ J can be 
achieved [l^[l3, H^l- In a typical experiment, atoms are 
cooled in an optical lattice of retroreflected diode lasers 
and then transferred into a magnetic trap where further 
cooling is to take place. This creates an arrangement 
of atoms where the resulting optical potential depths 
Vx^y^z are proportional to the laser intensities and can 
be expressed in terms of the recoil energy Efj = 
with m the atomic mass and k the wavelength num- 
ber. To prepare ID arrays two lattice lasers are given 
high intensities in such a way that hopping can only 
efficiently take place across one axis [loj. In terms of 
experimental parameters, hopping and repulsion coeffi- 
cients are given by J = A (^^^ g-C^yoZ-BR^;^ ^-^^ 



depth, V± the depth of the lattice in the transverse di- 
rections. Us the s-wave scattering length, d the lattice 
spacing and A, B and C fixed constants 0, [IJ]. Spatial 
variations in U and J can also be implemented using de- 
tuned lasers sent through specific sections of the lattice 
as reported in Ref. . 

From the Heisenberg equations of motion we know 



the general case non-linearities induced mostly by repul- 
sive terms keep us from getting reliable expressions for 
the ctk in terms of time. One fortunate instance is the 
repulsionless case where all T4 go to zero and we are left 
with 
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One can of course assume that all the coupling con- 
stants are equal, which corresponds to the usually stud- 
ied, and perhaps the most natural, setting of the Bose- 
Hubbard model. To contrast with this, one can also 
study a setting with non-uniform couplings to investi- 
gate whether better end to end entanglement can be ob- 
tained by appropriately engineering the couplings. As an 
example we will choose a non-uniform hopping distribu- 
tion that matches an angular momentum representation 
of length j = namely Jk = |V/c(iV- fc) (A = 2 for 

all the simulations presented in this work). The above 
couplings, when present in spin chains, are known to fa- 
cilitate a perfect quantum state transfer [2^, which has 
been developed in context of the idea of using spin chains 
to convey quantum states ^26*1. In a Hubbard model, 
the above coupling profile will transfer particles perfectly 
from one end of the chain to the other [12] (however, the 
aim of this paper is different, namely generating entan- 
glement between the ends). Moreover, because we are 
interested in studying end-to-end entanglement (EEE), 
hopping distributions are known to play a crucial role 
as it is particle tunnelling what determines how distant 
places get correlated with each other. So, in this work 
we compare results from two different hopping profiles, 
namely, the well known constant hopping (CH) Jk = 1 
and the above introduced perfect transmission hopping 
(PTH). The purpose in bringing in PTH is to explore 
what kind of physics is shown by a hopping profile with 
perfect transmission properties in situations where trans- 
port is either absent or not directly involved. The natural 
question in this context is whether perfect transmission 
properties, which involve dynamical synchronization at 
long scales, are in any way linked to the onset of quan- 
tum correlations among distant sites of a boson chain. 

On the other hand, when repulsive terms in Eq. ([1]) 
are considered we apply numerical methods. Here we 
make use of density-matrix techniques as presented in , 
also known as TEBD, combined with conservation and 
symmetry properties associated with the Hamiltonian. 
Briefiy, TEBD consists in writing the state in terms of 
canonical coefficients that are closely related to Schmidt 
vectors. These coefficients are updated through succes- 
sive applications of semi-local unitary operations that 
correspondingly only modify coefficients associated with 
the sub-spaces on which they act. This fact is then skil- 
fully exploited to work out an efficient way of simulating 
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both real and complex evolution. In addition to the stan- 
dard set of coefheients F and A employed in the canoni- 
cal representation of the state, we also retain the number 
of particles associated with every Schmidt vector in the 
decomposition. Memory storage is handled by making 
use of allocatable types and pointers in FORTRAN 95, 
which suits the mechanics of TEBD where the dimen- 
sion of the Hilbert space is updated adaptively. Usually, 
the TEBD coefficients that determine the local proper- 
ties of the chain are presented with explicit reference to 
every local dimension, i.e. r|^^', where « = 1, being 
M the number of total number of bosons in the chain. 
However, the only relevant component i above is com- 
pletely determined by the number of particles TOq and 
mp associated with Schmidth vectors \a) and in the 

expression ji/)) = Z]q./3,i r[^*|a) |i) through the rela- 
tion rrioi + mp ~\- i — M . Consequently, truncating the 
local Hilbert space as a way of attenuating the mem- 
ory requirement of the simulation is not necessary. Sim- 
ilarly, we use dynamical allocation and deallocation to 
store the canonical representation, which allows an ef- 
ficient handling of computer memory. The maximum 
number of Schmidt coefficients x as well as the size of 
other elements relevant to the simulation are updated at 
run time. After every updating step we retain all the 
Schmidt coefficients greater than A^^ x 10^^'', where A^^ 
is the greatest coefficient at site i. In ground state sim- 
ulations, X saturates and remains quite fixed until the 
state converges (Insets in Fig. |S]). For real time sim- 
ulations, on the other hand, this convenient saturation 
does not manifest as strongly as in the ground state case 
and in some cases, specially for long simulations, trunca- 
tion in the number of coefficients is important as a mean 
of sustaining the simulation efficiency. In our program 
we use second order Trotter expansion with time steps 
5t = where Nf. is the effective length of the chain, 

which is to be taken just as the number of sites in the 
chain unless something else be specified. Here we limit 
ourselves to Hamiltonians with symmetric distributions 
of parameters with respect to the centre of the chain and 
with an even number of sites and bosons. Additionally, 
all the simulations presented in this work correspond to 
M = N . Therefore, the system can be genuinely sim- 
ulated by focusing on the canonical coefficients of half 
of the chain plus the interaction between complementary 
half-chain blocks. The algorithm is applied recursively to 

simulate |?/'g) — UniT^ooTr — , until convergence to 

||e ^"\ipo)\\ 

the ground state \ipa) is accomplished according to the 
criterion |1 — (^(t)|V'(t + 6t))\ < 10~^^. Ground state 
simulations in repulsionless chains agree well with the 
exact ground state of Hamiltonian which for PTH 
chains reads 

with energy Eg — —XjN. Similar expressions for the 
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FIG. 2: Entanglement between both halves of the chain (main 
plots) and entanglement between the ends and the rest of the 
system (Insets) . In PTH chains, when repulsion is sufficiently 
strong tunnelling starts to take place throughout the whole 
chain, enhancing the amount of entanglement in the system. 
CH chains, on the other hand, display an entanglement reg- 
ister that can be understood as being originated from short 
scope hopping. 

ground state of CH chains can be found elsewhere. Ad- 
ditionally, real-time simulations coincide with available 
theoretical results [1, [H, H^l . 

As it was stated from the introduction of TEBD and 
further verified in several numerical studies [28,] , efficient 
simulations can be carried out as long as the state can be 
represented using reasonable low x, i-e., correlations are 
not too strong. As a reference one could take the infin- 
ity Hubbard model with unit filling, where such regime 
would correspond primarily to the Mott insulator phase 
where according to mean field theory ^ > 5.8 with 
z = 2d the number of nearest neighbours [3, UB] • An- 
other efficient simulation scenario corresponds to the case 
when, although hopping constants are strong and super- 
fluid features dominate, the Hamiltonian structure re- 
mains close to integrability, e.g. y ^ 1, since usually 
integrable dynamics takes place in subspaces. When the 
case is neither of the above mentioned, as it is for most of 
the simulations presented from now on, one can still see 
that X often remains well below the maximum allowed in 
the Hilbert space. 



III. RESULTS 

A. Ground state 

In order to calculate the entanglement between the 
ends of the chain, we need to get the corresponding re- 
duced density matrix. This task is by itself challenging, 
as the quantum state is given in terms of canonical co- 
efficients associated with TEBD and not in the standard 
basis. Similarly, the computation of such density matrix 
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FIG. 3: The emergence of EEE in the chain as the on-site 
repulsion or equivalently the fractional population at the end 
sites C, is increased. Linear behaviour can be observed inde- 
pendently of the chain size. Also, EEE in chains with PTH 
lingers on when the size of the chain is increased, in contrast 
to EEE in CH chains. 



is greatly improved in terms of speed and memory when 
exploiting number conservation. Once this density ma- 
trix is obtained, we have to quantify the entanglement 
using a measure which is appropriate for mixed states 
of two arbitrary dimensional quantum systems. So we 
choose the logarithmic negativity as defined in Ref. f2^. 
Moreover, to measure the entanglement between comple- 
mentary subsystems (30| , we use von Neumman entropy 
S = —tr[pAlog2{pA)) where pA represents the reduced 
density matrix of subsystem A. In this work we focus on 
chains with unit filling and open boundary conditions. 
Here we show results against the fraction of particles 

present on the ends C, = which allows a conve- 

nient depiction of the system phenomenology for differ- 
ent chain sizes. C varies according to the intensity of the 
repulsion in intermediate sites. When repulsion is small, 
few particles remain on the ends and ^ « 0. When re- 
pulsion is strong the amount of particles on the ends is 
maximum and C = | + • 

Irrespective of the hopping profile the ground state of 
repulsionless chains is highly superfluid, although with 
most of the hopping taking place in the central part, 
leaving the terminal positions nearly unoccupied. In this 
repulsionless regime entanglement is strong, but confined 
around the centre. For example, in these circumstances, 
it is strong between the left and right halves of the chain. 
When repulsion in intermediate positions is increased, 
particles are forced to hop through longer distances and 
correlations develop at longer scales. Similarly, bosons 
accumulate on the terminals according to the strength of 
the interaction in intermediate places. Notice that on ac- 
count of the repulsionless ends, particles can always find 
channels to hop, no matter how intense the repulsion in 
intermediate sites. Small repulsion leads to a decrease of 
entanglement between both halves and an increase of en- 



tanglement between both ends and the rest of the system, 
an indication that entanglement is being spread along the 
chain following the particle distribution profile. When 
repulsion is increased even more, the terminals start get- 
ting macroscopically occupied while the average number 
of bosons in intermediate sites falls asymptotically to- 
wards i , which means that strong tunnelling renders the 
quantum state of intermediate places into a state that 
involves superpositions of kets corresponding to one or 
zero bosons. In this case correlations among places near 
the ends and in opposite sides of the chain are strongly 
enhanced, in contrast to the correlations in the centre of 
the chain. This can be seen in the inset plots of Fig. [2] 
where von Neumman entropy between both ends and the 
rest of the system slightly comes down after the original 
redistribution of entanglement mentioned above occurs. 
The emergence of EEE in the chain is shown in Fig. [3l 
which shows that as repulsion in the centre becomes in- 
tense, the terminals will tend to entangle with each other 
more strongly than with intermediate sites. 

Once the conditions for the emergence of EEE are 
given, what determines how much entanglement can be 
generated is the hopping scope. Indeed, what is being re- 
ally adjusted when repulsion is turned on is the effective 
travelling length of bosons along the chain. As repulsion 
in the middle is tuned, bosons get to hop through longer 
distances and thereby their become more and more delo- 
calized. The behaviour displayed by von Neumman en- 
tropy between both halves in consistent with a regime in 
which entanglement is being continuously redistributed 
across the system as a result of increasing hopping scope, 
but the fact that there is no turning point in Fig. [2] 
for CH indicates that such hopping never takes place 
across the complete length of the chain. Nevertheless, 
such hopping is enough to induce EEE at finite iV, but 
not in the thermodynamic limit since EEE dies down 
against increasing chain size. In PTH chains, on the 
other hand, there exist a point in which hopping scope 
actually embraces the whole chain, which facilitates the 
particle transport from one end to the other, making the 
proportion of particles being held on the terminals use- 
ful for entanglement. Once this long scale hopping takes 
over, increasing the repulsion in intermediate sites rein- 
force end-to-end exchange and induces an increase in the 
amount of entanglement contained in the system, as can 
be seen in Figs. [5] and [31 At the same time, EEE be- 
gins to show a linear dependence with C, with a slope 
that is independent of the system size. In this case EEE 
remains even when the size of the chain is augmented, 
which means that long range correlations will manifest 
themselves in the thermodynamic limit. 



B. Dynamics 

Here we investigate how EEE entanglement shows up 
as a result of dynamics in a chain initially prepared in 
a Mott insulating state. The typical behaviour of entan- 
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FIG. 4: Entanglement between complementary blocks in a 
chain with A'' — 100 and U/ J = 0.4. In this simulation x = 50 
and 6t = 5 X 10"''. Entanglement saturates uniformly except 
close to the ends, where it displays smaller saturation values. 
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FIG. 5: Average number of bosons in chains with = 20 
and M = 20. x = 50. U = 100 for PTH (upper plot) and 
U = 200 for CH (lower plot). Stronger repulsion constants 
are necessary in PTH chains in order to maintain intermediate 
sites fermionized so that no more than one particle can be at 
these sites. The initial state is a Mott insulator with one 
boson per site. At the beginning bosons start to migrate 
towards the ends, leaving the intermediate of the chain with 
half a boson per site on average. As a result of the hopping, 
EEE emerges after some time (Fig[6ll. For longer times slow 
oscillations can be seen, but the ends remain well populated. 



glement in chains with homogeneous repulsion and hop- 
ping constants is shown in Fig. [H As hopping is bigger 
than repulsion, particle tunnelling generates rich dynam- 
ics and entanglement behaves non-trivially in contrast to 
the high repulsion regime. Entanglement saturation is 
identical for all sites except those very close to the ter- 
minals, where saturation takes place at smaller values. 
Such boundary effects are typical of open chains [1, |3l| . 
As a consequence, the ends are poorly entangled with 
the rest of system and between them as well. On the 



FIG. 6: EEE as a function of time for different chain sizes. 
X = 50. U = 100 for CH and U = 200 for PTH. The purpose 
is to restrain particle accumulation anywhere but the terminal 
sites. Entanglement in PTH chains is not only bigger but it 
also arises at early times than in CH chains. In both plots 
Log-negativity decreases with the size of the chain. Also, 
entanglement emerges well after bosons migrate to the ends 
(FigEJ. 



other hand, dynamics in repulsionless chains is known to 
lead to a progressive thermalization of reduced density 
matrices of each site [s^l and therefore any sort of two 
site entanglement is weak. In this sense, dynamic EEE 
in repulsionless chains is characteristically similar to that 
found previously in the ground state. As an alternative 
approach, we try the same methodology already applied 
and set all intermediate repulsion constants to high val- 
ues in order to artificially create a particle distribution 
favourable to EEE. The effect on the average number of 
particles can be seen in Fig. \5\ The ends get macro- 
scopically occupied and tunnelling in intermediate places 
intensifies. EEE arises some time after the ends have 
been occupied (Fig. [6]). In general, the bigger the chain 
the longer it takes for EEE to emerge. For PTH, a natu- 
ral time scale is determined by the transmission period, 
namely T — t: given our particular choice of A. Entan- 
glement arises after half of such period, just after infor- 
mation has got to travel from one end of the chain to the 
other, but there is not enough evidence that this is always 
the case, particularly for very long chains. Additionally, 
it is possible that entanglement will build up further for 
longer times, independently of the hopping profile, but 
this effect is difficult to observe as long time simulations 
require additional computational resources [H, [s^ . 



C. Perturbation approach and detection 

We would now ideally like to generate such an entangle- 
ment between the ends where the logarithmic negativity 
exceeds unity (the maximum value in a two qubit state) 
so that the possibilities of multiple occupation numbers of 



6 



60 

z 

60 

o 




FIG. 7: EEE in a chain with A'' = 14 and different values of 
A for both CH and PTH. PTH lines are thicker than their 
CH equivalents. EEE in the ground state can be improved as 
a result of dynamics in PTH chains, x = 200 for PTH and 
X = 100 for CH. 



where \'ijjg) is the ground state of if , A is a small num- 
ber and hj represents a local operator acting on site 
j. For consistency, we require [H ^^j^ihj] ^ 0. Here 
we consider that the local perturbations can be written 
as functions of the corresponding local number opera- 
tor, i.e. hj{Nj), where Nj = ajoj. We can establish 
a non-vanishing commutator between the perturbation 
and the Hamiltonian by choosing, for every site but the 
ends, hj{x) = kx, where k is the integer distance between 
site j and the closest end. Notice that this is equivalent 
to considering a chain with spatially dependent chemical 
potential in intermediate positions. Moreover, the pri- 
mary factor determining hi^^ = h is the optimization of 
boson-transfer from intermediate sites into the ends when 
an initial state Itpg) is set. Supposing (H) remains stable 
during perturbation evolution and correlations between 
the Hamiltonian and the perturbation can be ignored, 
we can use the number-of-boson basis. In this way, per- 
forming average-energy balances among sequential par- 
ticle distributions underlying an optimal trajectory we 
get 



ca 
< 




FIG. 8: Entanglement detection after bosons are sent through 
a 50:50 beamsplitter (see Ref. HJ]). The original state is 
extracted from the ground state of chains with N=20. Insets: 
Ground state convergence using the TEBD algorithm for PTH 
(left top) and CH (right bottom) in chains with A'' — 20. 



the ending sites is utilized highly beneficially. Looking at 
potential alternatives to the schemes discussed before, we 
now assume that a chain initially prepared in the ground 
state of a Hamiltonian with high intermediate repulsion 
and highly fermionized intermediate positions is set to 
evolve under the action of a small perturbation poten- 
tial corresponding to a physical mechanism that inter- 
nally pumps bosons in and out the ends. The idea is to 
use the perfect transmission properties at a perturbation 
level to generate entanglement as a result of transport. 
The quantum state at any time can be written as 



t{H+A , h, 



Hno + + h{no + -{- - k)) = — , (5) 

where no is the average number of bosons allocated in 
one end at t = 0. After some algebraic manipulations 
we find h{x) — C2x'^ + cix with ci = ( + 2rto)A and 
C2 = —A. With this perturbation, the mean number of 
bosons does not fluctuate much since the dynamics is still 
governed by a Hamiltonian with strong repulsion coeffi- 
cients in intermediate sites and bosons remain squeezed 
into the ends. The leading mechanism in the evolution is 
the equitable exchange of particles between the terminals 
and the rest of the system. In addition, one can expect 
that such exchange will improve the communication be- 
tween distant places and also the entanglement. Fig. [7] 
shows log-negativity as a function of time for different 
perturbation intensities. CH results are shown for com- 
parative purposes. In the initial stages of evolution the 
dynamics is characterized by an increase of quantum fluc- 
tuations on the ends of the chain accompanied by little 
change in the average number of bosons. Hence, EEE is 
enhanced while avoiding massive migration of bosons to- 
wards the centre of the chain. Significantly, entanglement 
generation is improved not by adding interaction terms 
to the Hamiltonian but by adding a local perturbative 
potential. 

Finally, we would like to comment about how this 
entanglement can be verified experimentally. Once the 
atoms condense in the ground state or after the dynam- 
ics has taken place, the detection scheme presented in [34 1 
could be used. Following such a scheme, bosons on the 
ends are sent through a 50:50 beamsplitter and then the 
number of particles on the outputs is counted. Entangle- 
ment can then be detected using eab = tr{alacp) ~ ^ ~ 
(4) tr{d\aNp) where tr(a\acp) is the number of particles in 
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one output and p is the reduced density matrix describ- 
ing the ends, cab = for separable states and cab > 
for entangled states. In this way, entanglement is char- 
acterized in terms of experimentally measurable param- 
eters. For illustrative purposes, we present in Fig. [8] 
a plot of eab obtained from the ground state of chains 
with PTH and CH. We conclude that cab correctly de- 
termines whether the state is entangled. In order to illus- 
trate how our simulations converge to the ground state 
and the computational cost involved we also include plots 
of X against the number of steps. Here we emphasise that 
in our program x defines the dimension of an allocatablc 
array which is being continuously resized according to 
the number of Schmidt coefficients in the canonical rep- 
resentation. 



IV. SUMMARY 

We have studied the ground state entanglement as 
well as the dynamical entanglement displayed by Bose- 
Hubbard chains for several configurations of parameters. 



Our results indicate that chains with strong repulsion 
coefficients in intermediate places are suitable scenarios 
for the emergence of long scale hopping that can lead 
to a subsequent development in long-range quantum cor- 
relations. Entanglement also emerges as a result of dy- 
namics in chains initially prepared in a Mott insulator 
state. In addition, we showed that ground state entangle- 
ment can be improved through perturbation dynamics. 
Irrespective of the issue of entanglement, the behaviours 
displayed by the ground and dynamical states that we 
have predicted for different repulsion profiles should be 
interesting to verify in an experiment. It would also be 
interesting to investigate how the transition from short 
scope to long scope hopping occurs. It could happen, 
and our results provide some evidence of this, that this 
transition is in fact a second order phase transition where 
EEE works as the order parameter. 
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